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Abstract 

The possibihty of quantum oscillations of the tunnel conductance and magnetoresistance induced 
by spin-wave excitations in a ferromagnet-ferromagnet-ferromagnet double barrier tunnel junction, 
when the magnetizations of the two side ferromagnets are aligned antiparallel to that of the mid- 
dle ferromagnet, is investigated in a self-consistent manner by means of Keldysh nonequilibrium 
Green function method. It has been found that owing to the s-d exchange interactions between 
conduction electrons and the spin density induced by spin accumulation in the middle ferromagnet, 
the differential conductance and the TMR indeed oscillate with the increase of bias voltage, being 
consistent with the phenomenon that is observed recently in experiments. The effects of magnon 
modes, the energy levels of electrons as well as the molecular field in the central ferromagnet on 
the oscillatory transport property of the system are also discussed. 
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I. INTRODUCTION 



In past decades, the spin-dependent transport properties in magnetic tunnel junctions 
(MTJs) have been extensively investigated both experimentally and theoretically, where a 



tigated botn experimentc 
.g. Refs. Q, 13,0, 0,3 for 



reviews) . It has been unveiled 



great progress has been made (see, e 
that owing to the conduction electron scatterings, the tunnel current through the MTJ is 
modulated by the relative orientation of magnetizations, giving rise to the so-called tunnel 
magnetoresistance (TMR) effect. As the quality of tunnel junctions is being improved, a 
large TMR, which is expected by practical applications, has been achieved in several systems. 
On the other hand, a reverse effect of TMR, coined as the spin transfer effect js, 7], has also 
been proposed, which predicts that the orientation of magnetization of free ferromagnetic 
layer can be switched by passing a spin-polarized electrical current, and spin waves could 
also be excited. This latter effect has been confirmed experimentally in a number of systems. 

Although single barrier MTJs already show abundant characteristics concerning the spin- 
dependent electrical transport, a double barrier magnetic tunnel junction (DBMTJ), in 
which the formation of quantum well states and the resonant tunneling phenomenon are 



theoretica 



13 



14 



15 



ly anticipated, has also attracted much attention in recent years 
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In order to observe the coherent 



tunneling thru the DBMTJ, people have attempted to improve the junction quality to elim- 
inate the influences from the interface roughness and impurity scattering, and remarkable 
advances have been achieved on this aspect. 

Recently, an unusual magnetotransport phenomenon in the ferromagnet-ferromagnet- 
ferromagnet (FM-FM-FM) DBMTJs was reported by Zeng et al. 27|]. They observed that, 
when the magnetization of center (free) magnetic layer was antiparallel (AP) to the mag- 
netization of the two outer (pinned) magnetic layers, the conductance and TMR oscillate 
distinctly with the applied bias voltage, while for the parallel (P) situation, no such os- 
cillation was seen. Unlike the previous oscillatory tunnel magnetoresistance, this unusual 
phenomenon can neither be explained by Coulomb blockade effect since the middle FM layer 
is continuous, the charge effect should be equal in P and AP configurations and the charging 
energy is negligibly small, nor be attributed to the resonant tunneling, because the observed 
period of oscillation is too small to account for the energy level spacing of the quantum well 
states. Considering that the conductance oscillation is asymmetrical for P and AP configu- 
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rations, and the energy level of the unusual phenomenon is the same as the typical energy of 
a magnon, one may speculate that the unusual oscillation behavior could be induced by the 
magnon-assisted tunneling 23]. This is because in the AP state, the nonequilibrium spin 
density, which is proportional to the applied bias, could be accumulated near the interfaces 
in the middle region to emit spin waves, and the magnon-assisted tunneling would contribute 
to the conductance, while in the P state the spin wave emission is forbidden due to the spin 
angular momentum conservation, as discussed previously 



As there is no previous theoretical study devoting to the investigation on the possible 
quantum oscillations induced by spin wave excitations, in this paper, by using the nonequi- 
librium Green function method, we shall examine theoretically the above-mentioned idea 
by studying the possibility of magnon-assisted tunneling in the FM-FM-FM DBMTJ, and 
explore whether the magnon-assisted tunneling could really cause the oscillations of the 
differential conductance and TMR with the applied bias voltage. 

The rest of this paper is organized as follows: In Sec. II, a model is proposed. The 
tunnel current and relevant Green functions are obtained in terms of the nonequilibrium 
Green function technique in Sec. III. In Sec. IV, the transport properties of the system are 
numerically investigated, and some discussions are presented. Finally, a brief summary is 
given in Sec. V. 

II. MODEL 

Let us consider a FM-FM-FM DBMTJ with three FM layers separated by two thin 
insulating films. Suppose that the left (L) and right (R) FM electrodes with magnetizations 
aligned parallel are applied by bias voltages —Vjl and V jl^ respectively. The magnetization 
of the middle FM layer is presumed to be antiparallel to those of the L and R electrodes so 
that spin waves can be emitted in the middle FM layer because of spin accumulation. The 
schematic layout of this system is depicted in the inset of Fig. 1(a). The Hamiltonian of 
the system reads 

H = Hl + Hr + Hc + Hlc + Hcr, (1) 
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with 

Ha = Y^ Sk^aat^^ak^a, {a = L, R) (2) 



kacr 



He = ^£ka4„Cka + "^fhJgb+bg, (3) 
ka q 



Hlc = J2^k,kK.c,^ + h.c.) 

ki^ka 

kLkq 



kLkq 



HcR = J2Kk«.^k. + h.c.) 

kjika 



kukq 



)^^Ti^kq^{<4R^Ckiii^ + CfcT^fcfli^g + (^Ri^kA + 4|afc^t6g), (5) 



where ak^a and Cfco- are annihilation operators of electrons with momentum k and spin a in 
the a electrode and in the middle FM layer, respectively, Sk^a- = £kc, ~ crMa — eVa with Sk^ 
the single-electron energy and the molecular field in the a electrode, Ska = ^k — crM 
with Ek the single-electron energy and M the molecular field in the middle FM layer, bq 
is the annihilation operator of magnon with momentum q in the middle region, hujq is the 

magnon energy, = with = b^bq is the number of magnons, S{q) = 3 — 11^ where 

q 

S" = 1/2 is the spin of electron, T^^^ are tunneling matrix elements of electrons between the 
a electrode and middle FM layer, T^^^^ are coupling matrix elements between the electrons 
in a electrode and magnons in the middle FM region. 

It is noting that H^c (Hcr) describes the coupling between electrons in the L (R) elec- 



trode and electrons as well as magnons in the central FM region, w. 
"^kakq (H]) and ([5]) are due to the s-d exchange interactions 

erahty, we further assume T^^^ = 7T^^^ in the following discussions. 



lere the terms containing 



28|. Without loss of gen- 
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III. TUNNEL CURRENT AND GREEN FUNCTIONS 



A. Tunnel Current 

Starting from Eq. ([T]), after some cubersome but straightforward calculations, one may 
obtain the tunnel electrical current / 

I = h^ + h,, (6) 



kf^k kj^kq 



I2S 



'J^fy^Jkr.kq^'kl.kTA^^ t)]^ 



khkq 



kijk ki^kq 



'2S 



+ \/ 'T7y^JkT.kqGk^k,.\{^^ ^)]^ 



k^kq 



where the lesser Green functions are defined as 



Gtk,.M = M,At'yUt)), 



(7) 



G&,i(t,t')=^«i(0&.(OcfcTW)- (9) 

It should be remarked that in the above derivations, we have made decoupling approxima- 
tions for the terms containing to simplify the calculations. From these above equations, 
one may see that to get the tunnel electrical current, the lesser Green functions must be ob- 
tained. In the following, we shall employ Keldysh's nonequilibrium Green function method 
to get all self-consistent equations to determine the lesser Green functions. As the lesser 
Green function is closely related to the retarded and advanced Green functions according 
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to Keldysh formalism, the relevant retarded and advanced Green functions of electrons and 
magnons should be first calculated. 

Accordingly, the differential tunnel conductance (G) is obtained by G{V) = dI{V)/dV, 
and the TMR can be calculated by TMR = (1 — Gij/Gn) where G-^ is the differen- 

tial conductance when the magnetizations of the middle FM and the side FM are aligned 



antiparallel (parallel) . 

B. Green Functions of Electrons 



Let us define useful retarded Green functions for electrons as 

GLk'aM = -^0{t-mCka{t),4,,{t')}), (10) 

GLk.At,t') = -iOit-t'){Mt),al,it')}), (11) 

GI^LM = -i^(i-t')(KTW,<i(^')&.(^')}), (12) 

Gl^'liM = -^^(^-0({c.TW,4i(0&.(^')}), (13) 

GtLiM = -te{t-t'){Mt),al,{t')b^m, (14) 

Gl^Ut^t') = -ie{t-t'){Mt),4,mm. (15) 



In terms of the equation of motion, after a tedious calculation, up to the third-order of 
Green functions, we get the following equations 

- ^k^dGh,^^{s) = r + ^E('^-K))]E^^^m(^) 

+/f^'^E^S.(^)' 

k'q 
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" a ka 



N 



a kaq 



l2S 



I q' k" ) a ka 
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k" 



^2^ 



E(^-K))]EE^^i^a(^) 



a feo 



'25 



'25 



(£-£,,,-K)G^.i.'T(^) = fr'+^E(^-K')) + E(4TC^"T)] EE^t 

I ^ q' k" ) a kc 



'25 



+V aT^'^^^^^ + ^ " E(4'tC^"t))EE^^^t.«t(^)' 



k" 



a ka 
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where we have presumed, for simphcity, the couphng matrix elements T^^^^, and T^^^ in- 
dependent of momentum by considering that only those electrons near the Fermi surface 
participate in the transport process, and n^^o- = al^^a^^^. 

From these equations, the required Green functions can be obtained self-consistently. 
On the other hand, the lesser self-energy can be approximated by Ng's ansatz 29|: 
S< = S<(S5 - - S"^), where - = G^-^ - G"-^, SJ, and S< are given by the 

following equations 



-0TTV^7 ^oin^' I ^ j 2 ~ ^ 



(16) 



So<^(5) So<^(5) \f ^TL^fie-f)+^Tn^f{e+f) \ 

S<^(5) S<^(5) y 1^ ^rLjie-f)+tTnJie+f) ) 

where Taa{^) is the linewidth function defined by Taa{^) = 27r^po-(A;Q,) j?"^^^!^ with Pa{ka) 
the density of states of electrons with momentum ka and spin cr in the ath FM electrode, 
and f{e) is the Fermi distribution function. By means of = CS^G", the lesser Green 
functions can be procured. 



C. Green Functions of Magnons 

As the number of magnons, A^, enters into the formalism, we need to obtain the Green 
functions of magnons to determine self-consistently. Define the retarded Green function 
of magnons as 

G;,(t,t') = -^^(t-t')([6,(t),6+(t')]). (17) 
By using the equation of motion, we have 



A^ 



^ E T^TLa-G^^^) + Gt\^) - Gt\e) + Glfie)], (18) 
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where G^j^^^'(s) (i — 1, • • • ,5) are the Fourier transforms of the Green functions defined as 
below 



GZ\t,t') = -ie{t-t'){[b,(t),atJt'Mt')]), 

Glfm = -^^(^-^')([&.(t),&+(t')ati(^>^i(^')])' 
Glf{t,t') = -ie{t-t'){[b,{t),b+{t)ct,{t)aUt)])- 



By using repeatedly the equation of motion, and making appropriate cut-off approximations, 
up to the second order, we have 



e — hWq — Ek^l + Ski 

(21) 



e — hlOq + Ek^i — £fe-f 

(22) 



e — hxjg — Ek^i + Ski 

(23) 



(24) 




(19) 




(20) 
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The number of magnons can thus be obtained by the spectral theorem 

^ = EK) = / = E / - Gl^{e)l (25) 

9 9 9 

where /^(e) is the Bose distribution function. 

To get physical quantities under interest, all above equations should be numerically solved 
in a self-consistent manner. 



IV. RESULTS AND DISCUSSIONS 



To proceed the numerical calculations, we need to make some assumptions. Since the 
number of the above self-consistent equations nonlinearly increases with increasing the num- 
ber of wave vectors of electrons and the number of spin-wave modes in the middle FM 
region, which makes the calculations too complicated to perform, for the sake of simplic- 
ity but without losing the generality, we shall only consider the situations where both the 
numbers of k and q taken in the following calculations are not so large that the numeri- 
cal calculations can be readily proceeded. This is plausible, because the magnon-assisted 
transport property mainly depends on the low-lying quantum well states of electrons in the 
middle FM, and only the lower modes of spin waves are easy to emitj^, leading to small 



energy levels of magnons 



271]. Besides, considering that only those electrons near the Fermi 



surface participate in the tunneling process, we may take Ska ~ ^kp — o'M, denoted by 
and for spin up and down electrons, respectively. In addition, we suppose that the 
two side FM electrodes are made of the same materials, i.e.. Ml = Mn, Pl = Pr = P, 
where Pl(r) = [Tl{r)^ - ^LiR)i]/[^L{R)^ + ^l{r)i] is the polarization of the left (right) FM 
layer. Then, the linewidth function can be written as T^-^^^ = F^^^ j = ro(l ± P), where 
Fo = Fj^(/j)|(P = 0) = Fi(/j)j(P = 0) will be taken as an energy scale. In the following, we 
will take P = 0.7, ksT = 0.04Fo, Iq = ^ and Gq = ^ will be taken as scales for the tunnel 
current and the differential conductance, respectively. 



A. Effect of Magnon- Assisted Tunneling 

In order to study whether the quantum oscillations of the conductance and TMR ob- 
served in the FM-FM-FM tunnel junction are induced by spin-wave excitations owing to 
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spin accumulation, let us first examine the bias dependence of the transport properties by 
considering the effect of magnon-assisted tunneling in the AP state. 

For given energy levels of the electrons and magnons in the middle FM region, the bias 
dependent tunnel current (///q), the differential conductance (G/Gq) and TMR for different 
7 (= T"^ /T'^) are shown in Fig. 1. It is observed that the tunnel current increases nonliniearly 
with increasing the bias voltage. The twisted behavior of I/Iq for 7 7^ 0, presented in Fig. 
1(a), comes from the magnon-assisted tunneling, as no such behaviors are found for 7 = in 
either AP or P state. This can be clearly seen from the bias dependence of the differential 
conductance, shown in Fig. 1(b), where the peaks and dips appear for appropriate 7. 
Correspondingly, the TMR shows oscillating behavior with increasing the bias voltage, as 
demonstrated in Fig. 1(c). This observation manifests that the quantum oscillations of 
the conductance and TMR in the FM-FM-FM tunnel junction can be caused by spin- wave 
excitations, because when we turn off the effect of spin-wave excitations, the oscillating 
behaviors of G and TMR disappear. Note that the peak in G/Gq and one dip and one 
peak in TMR are from the quantum resonant tunneling of electrons. A larger TMR can be 
obtained for large 7, and due to the s-d exchange interactions that could lead to spin-flip 
scatterings, the TMR can be negative, as presented in Fig. 1(c). 

It should be remarked that the above oscillating behaviors of G and TMR appear only 
when 7 is in a suitable range, say, when 7 is too small, no oscillations can be observed, while 
7 is too large, the self-consistent equations have no solutions. 

The reason for the appearance of oscillations is that, when the applied bias voltage 
exceeds a certain value, the non-equilibrium spin density can be accumulated in the middle 
FM region, and spin waves can be excited. When polarized electrons from the left FM layer 
tunnel into the central FM layer, they are subject to scatterings from not only the polarized 
electrons in the central region but also the spin waves from the accumulation owing to s-d 
exchange interactions. It is possible that the electrons may turn their spin directions by 
emission and absorption of magnons, leading to that the tunnel conductance and TMR 
oscillate under a combination of effects of magnon-assisted tunneling as well as quantum 



resonant tunneling 



23| through the quantum well states in the central FM region. 



Apparently, the polarized electrical current can excite magnons, while these magnons 
participate in the tunneling process and in turn influence the tunnel current. The number of 
magnons must be estimated self-consistently. As an example, in Fig. 2, the bias dependent of 
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2 3 4 5 6 7 8 



FIG. 1: (Color Online) Bias dependence of (a) the tunnel current I, (b) the differential conductance 
G, and (c) TMR for the different 7, where hiOg = 0.5To, = 3.0Tq and = 2.0ro. 

the number of magnons for some parameters is presented. We may see that the number of 
magnons oscillates with the bias voltage, which could be the main reason for the oscillatory 
transport property of the system. 



B. Effect of Magnon Modes 

There are a number of factors including magnon energy hwg that can affect the transport 
behavior of the FM-FM-FM tunnel junction. The bias dependence of the current, the 
different conductance and the TMR for different magnon energies is shown in Fig. 3. It can 
be found that with increasing huq, apart from some quantitative changes of peak positions 
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FIG. 2: (Color Online) Bias dependence of the number of magnons N in the AP state. 

and amplitudes, there are not much quahtative changes of the current, conductance and 
TMR. Therefore, for a given magnon mode, the magnon energy does not affect quahtatively 
the transport oscillating behavior of the system. It is noted that, while some positions of 
the peaks and dips of G and TMR are influenced by the magnon energy, the others are not. 
This observation indicates that the magnon energy is only one of factors determining the 
positions and amplitudes of the peaks. 

For the sake of simplicity, in the aforementioned analysis we have adopted a single spin- 
wave mode. Whether are the transport properties of the system much affected qualitatively 
when we take more spin-wave modes into account? The answer is presented in Fig. 4, 
where we have taken two and three spin-wave modes to get the tunnel current, differential 
conductance and TMR. For a comparison, we have also included the case with single mode. 
One may see that at low biases, the magnon modes do not have so much effect on the 
behaviors of ///q, G/Gq and TMR, but at higher voltages the magnitudes of the current, 
differential conductance as well as TMR change somewhat remarkably. This is because 
at low biases the spin accumulation effect is small, and the interaction between tunneling 
electrons and spin-wave modes is weak, leading to the transport properties less influenced; 
at large biases the spin accumulation effect becomes more pronounced, and the interactions 
between electrons and magnons are strong, the transport behaviors of the system are thus 
altered quantitatively. Note that the round peaks and dips in the curves of the tunnel 
current shown in Fig. 4 may come from combinations of the magnon-assisted tunneling as 
well as the quantum resonant tunneling. 
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FIG. 3: (Color Online) Bias dependence of (a) the tunnel current I, (b) the differential conductance 
G and (c) TMR for the different fiWq, where 7 = 0.2, = S.OFq and = 2.0ro. 

C. Effect of Electron Level in the Middle FM Layer 



The bias dependence of the current, the differential conductance and the TMR for dif- 
ferent energy levels {e^,^ = e) of the electrons in the middle FM region is given in Fig. 5. It 
can be observed that with lifting the energy levels of electrons in the central FM layer, the 
peak and dip positions of I/Iq, G/Gq and TMR change dramatically with increasing the 
bias voltage, while the shapes of the curves retain quite similar for different energy levels 
of electrons in the central region. This signifies that the energy levels of electrons in the 
middle FM region do not affect the oscillating behavior itself of the transport properties. 
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FIG. 4: (Color Online) Bias dependence of (a) the tunnel current I, (b) the differential conductance 
G, and (c) TMR for different spin- wave modes, where 7 = 0.2, e-f = S.OFq and ej = 2.0ro. 



but affect the positions of oscillating peaks and dips. This can be understandable, because 
the transport behavior of electrons are mainly determined by the scatterings from polarized 
electrons and magnons to that the electrons are subject, when the energy levels of electrons 
in the middle region are promoted, the resonant energies in magnon-assisted and quantum 
resonant tunneling processes become different, resulting in the behaviors shown in Fig. 5. 



D. Effect of Molecular Field in the Middle FM Layer 

The molecular fields of the middle FM layer could also have effect on the transport 
properties of the system. The bias dependence of the current, differential conductance 
and TMR for different molecular fields of the central FM region is shown in Fig. 6. It 
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FIG. 5: (Color Online) Bias dependence of (a) the tunnel current I, (b) the differential conductance 
G, and (c) TMR in the AP state for different energy levels of electrons in the middle FM region, 
where 7 = 0.2, fvjjq = O.SFo and M = LOFq. 

is seen that as the molecular field increases, the magnitude of the tunnel current becomes 
smaller, and the oscillations of the differential conductance as well as TMR become more 
apparent, where not only the number but also the positions of the oscillating peaks change 
with increasing the molecular fields. This fact suggests that the level spacing between the 
majorty and minority subbands of electrons in the central region plays an important role in 
the oscillating behaviors of the differential conductance and TMR in the present system. 
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FIG. 6: (Color Online) Bias dependence of (a) the tunnel current I, (b) the differential conductance 
G, and (c) TMR for different molecular fields of the middle FM layer, where 7 = 0.2, e = 2.5ro 
and hLOq = O.SFo. 

V. SUMMARY 

In summary, we have probed the possibility of quantum oscillations of the differential con- 
ductance and TMR in the FM-FM-FM tunnel junction in the AP state. By self-consistently 
taking the s-d exchange interactions between conduction electrons and the nonequilibrium 
spin density induced by spin accumulation in the middle FM layer into account, we have 
found that the differential conductance and TMR indeed oscillate with increasing the bias 
voltage, thereby theoretically confirming qualitatively the inferrer and experimental results 
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presented in Ref. 27|]. It has been unveiled that the average number of magnons oscillates 
with the bias, which could be the main reason for the oscillations of the conductance and 
TMR of the system. When we turn off the s-d exchange interactions, i.e., taking 7 = 0, 
no oscillations of the conductance and TMR were observed, showing that the oscilations 
are indeed caused by the spin-wave excitations induced by spin accumulations. In the P 
state, owing to the absence of spin accumulation 2l|, |22|, no oscillations of the conductance 



and TMR with the bias can be found. We have also investigated the effects of the magnon 
modes, the energy levels of electrons as well as the molecular field in the middle FM region, 
and found that in spite of changes of the positions and amplitudes of the oscillating peaks 
and dips, the oscillatory behavior of the transport properties is not qualitatively affected. 
We anticipate that our findings could offer clues for better understanding the experimental 
observation presented in Ref. {2?! . 

Finally, we would like to remark that our preceding discussions could be applicable to 
the system with a magnetic quantum dot coupled to two ferromagnetic electrodes, where 
the oscillatory behavior of the transport properties with the bias would be expected if the 
spin accumulation effect is not neglected. The work toward this direction is under progress. 
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